Absence of an Almeida-Thouless line in Three-Dimensional Spin Glasses 
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We present results of Monte Carlo simulations of the three-dimensional Edwards- Anderson Ising 
spin glass in the presence of a (random) field. A finite-size scaling analysis of the correlation length 
shows no indication of a transition, in contrast to the zero-field case. This suggests that there is no 
Almeida-Thouless line for short-range Ising spin glasses. 
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Since the work of Ballesteros et al. [l| , there has been 
little doubt that a finite-temperature transition occurs in 
three-dimensional spin glasses 2J. However, the behavior 
of a spin glass in a magnetic field is less well understood. 
In mean field theory |2(, which is taken to be the solu- 
tion of the infinite-range Sherrington-Kirkpatrick (SK) 
model 0, an Ising system has a line of transitions 
in a magnetic field known as the Almeida-Thouless 
(AT) line. This line separates the paramagnetic phase at 
high temperatures and fields from the spin-glass phase 
at lower temperatures and fields. Although there is no 
change of symmetry at this transition, the relaxation 
time diverges (and for short-range systems so does the 
correlation length as we shall see) . In the spin-glass phase 
below the AT line, there is "replica symmetry breaking" 
in which the free energy landscape breaks up into differ- 
ent regions separated by infinite barriers, and the distri- 
bution of relaxation times extends to infinity. 

It is important to know whether the AT line also occurs 
in more realistic short-range models, since the two main 
scenarios that have been proposed for the spin-glass state 
differ over this issue. In the "droplet picture" 0,01 El E3 
there is no AT line in any finite-dimensional spin glass. 
By contrast, the "replica symmetry breaking" (RSB) pic- 
ture 0,0,0,0] postulates that the behavior of short- 
range systems is quite similar to that of the infinite-range 
SK model which does have an AT line as just mentioned. 
Both scenarios are illustrated in Fig. ^ 

Experimentally, it has proved much more difficult to 
verify the transition in a field than for the zero- field tran- 
sition. For the latter, the divergence of the nonlinear sus- 
ceptibility gives clear experimental evidence of a transi- 
tion, but unfortunately this divergence no longer occurs 
in a magnetic field. Experiments have therefore looked 
for a divergent relaxation time, and a careful analysis by 
Mattsson et al. [15J finds that this does not occur in a 
field. However, not all experimental work has come to 
the same conclusion 

In simulations, it is most desirable to perform finite- 
size scaling (FSS) on dimensionless quantities for rea- 
sons that we will discuss below. One such quantity, the 
Binder ratio, gave some evidence for the zero-field tran- 
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FIG. I: (a) H—T phase diagram expected according to RSB 
for the short-range case. For T < T C (H) there is a spin- 
glass phase (SG), whereas for T > T C (H) the system is in 
the paramagnetic (PM) state. The value of the critical field 
at T = is called Hat- (b) H-T phase diagram following 
the predictions from the droplet picture. A spin-glass phase 
exists only for H = 0. 



sition 0,0]. However, the Binder ratio turns out to be 
very poorly behaved in a field |l9l | in short-range systems, 
while for the SK model it does indicate a possible transi- 
tion [2(J, although not with any great precision. Results 
of out of equilibrium simulations on large lattices in four 
dimensions 0] were interpreted as evidence for RSB, al- 
though it is not completely clear that the true equilibrium 
behavior is probed by this procedure 0]- By contrast, 
simulations [23j corresponding to experimental protocols 
in (non-equilibrium) aging experiments have been ana- 
lyzed in terms of a "dynamical crossover" consistent with 
the droplet picture. 

Houdayer and Martin 0] carried out interesting cal- 
culations at T = to determine Hat, the critical field at 
T = 0, see Fig. ^ f° r a simple cubic lattice in three 
dimensions. Their results indicated that Hat = 0, 
i.e., there is no AT line, although a subsequent zero- 
temperature study by Krzakala et al. 0] found some 
evidence of a critical field for H ~ 0.65, which is much 
less than the "mean field" value for this lattice (2|| of 
around 1.86. However, Krzakala et al. could not ex- 
clude the possibility that the critical field is zero. 

As noted above, the (dimensionless) Binder ratio has 
provided some evidence for zero-field transition at finitc- 
T. However, a much sharper signature of the zero-field 
transition in three dimensions is provided by the corre- 
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lation length Q from which a dimensionless quantity is 
formed by dividing by the system size L. This approach 
should also provide evidence from equilibrium calcula- 
tions for a transition in a field, if one occurs, and in 
this paper we use it to determine whether there is an 
AT line in a three-dimensional Ising spin glass. Our con- 
clusion will be that there is not, at least down to fields 
significantly smaller than the value of 0.65 suggested by 
Krzakala et al. |2fij . 

The Hamiltonian we study is given by 



(1) 



in which the Ising spins Si — ±1 lie on the sites of a 
simple cubic lattice of size N — L 3 (4 < L < 12) with 
periodic boundary conditions, and the nearest neighbor 
interactions Jy are independent random variables with a 
Gaussian distribution with mean zero and standard devi- 
ation unity. At each site there is a field hi which, like the 
bonds, is randomly drawn from a Gaussian distribution, 
and whose mean and standard deviation are given by 



[hiU = 0, [hf]\i 2 = H T , 



(2) 



where [• • -] av denotes an average over the disorder. For 
a symmetric distribution of bonds, the sign of hi can be 
"gauged away" so a uniform field is completely equivalent 
to a bimodal distribution of fields with hi — ±H. Our 
choice of a Gaussian distribution, which still has an AT 
line in mean-field theory, also puts disorder into the mag- 
nitude of the hi. We use a Gaussian distribution, rather 
than a uniform field, in order to apply a very helpful test 
for equilibration, discussed below. 

To determine the correlation length we calculate the 
wavevector-dependent spin-glass susceptibility which, for 
nonzero fields, is defined by 



ik-(Ri-Rj) 



(3) 

where (••j)t denotes a thermal average. As in earlier 
work |]J, |27j the correlation length of the finite system is 
defined to be 
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(27r/L,0,0) is the smallest nonzero 



where k m j n 
wavevector. 

Now £l satisfies the finite-size scaling form 



^- = x(l x i»[t-t c (h t )] 



(5) 



where v is the correlation length exponent and T c (H r ) is 
the transition temperature for a field strength H T . Note 



TABLE I: Parameters of the simulations for H r = 0.3. N sarap 
is the number of samples, N avreep is the total number of Monte 
Carlo sweeps for each of the 4Nt replicas for a single sample, 
r m i n is the lowest temperature simulated, and Nt is the num- 
ber of temperatures used in the parallel tempering method. 
For other values of H r , we used the same parameters but only 
simulated L = 4, 6, and 8. 
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that there is no power of L multiplying the scaling func- 
tion X, as there would be for a quantity with dimensions. 
This greatly simplifies the analysis since the critical point 
can be seen by inspection as the temperature where data 
for different sizes intersect. 

On the AT line, T = T C (H T ), the "connected correla- 
tion function" (SiSj)T — (Si)T(Sj)T becomes long range 
and so, for an infinite system, the correlation length and 
X SG (0) diverge while for a finite system, £l/L is inde- 
pendent of L. Below the AT line, according to RSB the 
correlation functions no longer have a "clustering prop- 



erty", 



lim 
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5 ((SiS,-)T - (SiWiSjh) + 0, so 



X SG (0) oc L d and £,l/L increases with L. Hence, ac- 
cording to RSB, the behavior of should be quali- 
tatively the same as at the zero-field transition, namely 
it decreases with increasing L above the transition, is 
independent of L at the transition, and increases with 
increasing L below the transition. 

We use parallel tempering to speed up the simulations 
but unfortunately it is less efficient in a field than in 
zero field |28|, [29( , because "chaos" with respect to a field 
is stronger than chaos with respect to temperature. As 
a result, the computer time increases very rapidly with 
increasing L, so it is unlikely that we will be able to study 
larger sizes in the near future without a better algorithm. 
In order to compute the products of up to four thermal 
averages in Eq. 10 without bias we simulate four copies 
(replicas) of the system with the same bonds and fields 
at each temperature. 

Parameters of the simulation are shown in Table [I] 
Most of our work is for Hr = 0.3 since this is smaller 
than the predicted [3, I^Ol zero temperature value of 
Hat = 0.65, but is not so small that the results would 
be seriously influenced by the zero-field transition. The 
lowest temperature is 0.23 well below the zero- field tran- 
sition temperature which is about [3JJ 0.95. 

For a Gaussian distributions of bonds and fields, the 
expression for the average energy, U = [(W)r]av, can be 
integrated by parts with respect to the disorder distribu- 
tion, with the result 



U = U(q,qi) 



(6) 
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FIG. 2: An equilibration plot for L = 6, H T = 0.3, T = 0.3 
showing that the data for the average energy U and the quan- 
tity U(q, qi), defined in Eq. Q, approach their common equi- 
librium value from opposite directions as the number of Monte 

Carlo SW66pS -/Vswecp 

increases. The inset shows data for £,l/L 
indicating that it has equilibrated when U and U(q,qi) have 
become equal. 




FIG. 3: Data for £ L /L for H T = for different sizes. Note 
that there are clear intersections at roughly a common tem- 
perature, with the data splaying out at lower temperatures. 
The temperature of the intersections is the zero-field transi- 
tion temperature, marked T c in Fig. Q In this data, 5000 
samples are used for the disorder average in each system size. 



where 



U(q,qi) = -— h -t^Hj. , 
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where z (= 6 here) is the number of neighbors, q is the 
spin overlap given by 



1 N 



and qi is the "link overlap" given by 
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zN 



T lav 
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In Eqs. (HJ and 10, "(1)" and "(2)" refer to two copies 
of the system with the same bonds and fields. Because U 
will decrease as the system approaches equilibrium and 
q and qi will increase (since we initialize the spins in 
the two copies in random configurations), U and U(q,qi) 
approach their common equilibrium value from opposite 
directions and so Eqs. © and Q can be used as an 
equilibration test. This is a g eneralization to finite fields 
of a test used previously [33 . Figure |21 shows that these 
expectation are born out. We accept a set of runs as 
being equilibrated if U = U(q,qi) within the error bars. 
The inset to the figure shows that has equilibrated 

when U and U{q,q{) have become equal. 
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FIG. 4: Data for £ L /L for H r = 0.3 for different sizes. Note 
that in contrast to the zero-field data in Fig. |!|]there is no sign 
of intersections down to the lowest temperature T = 0.23. 
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It is useful to compare results in a field with those at 
the zero-field transition. Hence in Fig.[3]we show data for 
£l/L f° r H T = for sizes up to L = 12. For these results 
we set (Si)T = in Eq. There are clear intersections, 
with data splaying out at lower temperatures, indicating 
a transition at T = T c = T c (H r = 0), in the region 0.95- 
1.00, in agreement with Marinari et al. |3l|. 

However, the analogous results for H r = 0.3 shown 
in Fig. 0] have no sign of an intersection for sizes up to 
L = 12 at temperatures down to T = 0.23, which is 
considerably below the zero-field transition temperature 
of about 0.95. This provides quite strong evidence that 
there is no AT line, except possibly for fields less than 0.3. 
In order to test this possibility we have also performed 
simulations down to H r = 0.05 (for 4 < L < 8), and again 
found no intersections. We have also performed simula- 
tions in a uniform field, finding that data are very similar 
to those for the random fields, and have no intersection 
down to the lowest field studied, H = 0.1. 

To go to low fields without passing too close to the 
zero-field transition, we also tried a diagonal "cut" in 
the H Y -T plane with H r /T kept fixed at the constant 
value of 0.7. However, the equilibration problems were 
even more severe than for H r fixed at 0.3, and so we have 
not been able to get useful data for this case. 

To conclude, our finite-temperature Monte Carlo sim- 
ulations provide simple, direct evidence from equilibrium 
calculations that there is no AT line in three dimensions. 
Of course, the numerical data cannot rule out a transi- 
tion at exceptionally small fields, or the possibility of a 
crossover at much larger sizes to different behavior, but 
we see no particular reason for these scenarios to occur. 
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